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(57) Abstract: The invention 
relates to a method, an X-ray 
apparatus and a computer 
software product, which make 
it possible to determine the 
distributions of material density 
and atomic number in objects 
to be examined of the human 
or animal type and in objects 
to be examined that are used in 
material and safety engineering. 
The functional dependency of 
a first X-ray absorption value 
of density and atomic number is 
determined in the instance of a 

Z first X-ray spectrum, and at least 
the functional dependency of a 
second X-ray absorption value of 

^ = density and atomic number is determined in the instance of a second X-ray spectrum. Based on a recording of a first distribution of 
X-ray absorption values of the object to be examined in the instance of a first X-ray spectrum, and on a recording of at least one 
second distribution of X-ray absorption values of the object to be examined in the instance of a second X-ray spectrum, the values 
for density and atomic number are determined by comparing the functional dependency of a first X-ray absorption value of the first 
distribution of X-ray absorption values with the functional dependency (ies) of the X-ray absorption values, which are assigned to 
the first X-ray absorption value, of the second and/or other distributions of X-ray absorption values. 




(57) Zusammenfassung: Beschrieben werden ein Verfahren, eine Rontgenapparatur und ein Computersoftwareprodukt, die ein 
Ermitteln der Verteilungen von Materialdichte und Ordnungszahl in Untersuchungsobjekten menschlicher wie tierischer Natur und 
in Untersuchungsobjekten der Material- und Sicherheitstechnik ermoglichen. Es wird die funktionale Abhangigkeit eines ersten 
Rontgenabsorptionswertes von Dichte und Ordnungszahl bei einem ersten Rontgenstrahlspektrum und zumindest die funktionale 
Abhangigkeit eines zweiten Rontgenabsorptionswertes von Dichte und Ordnungszahl bei einem zweiten Rontgenstrahlspektrum 
bestimmt. Ausgehend von einem Aufzeichnen einer ersten Verteilung von Rontgenabsorptionswerten des Untersuchungsobjekts bei 
einem ersten Rontgenstrahlspektrum und einem Aufzeichnen zumindest einer zweiten Verteilung von Rontgenabsorptionswerten des 
Untersuchungsobjekts bei einem zweiten Rontgenstrahlspektrum werden die Werte fiir Dichte und Ordnungszahl 
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aus einem Vergleich der funktionalen Abhangigkeit eines ersten Rontgenab sorptions wertes der ersten Verteilung von Rontgenabsorp- 
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Bes chr e i bung 

Verfahren zur Ermittluhg von Dichte- und Ordnungszahlvertei- 
lungen bei radiographischen Untersuchungsverfahren 

5 

Die vorliegende Erfindung betrifft ein Verfahren, eine Ront- 
genapparatur und ein Computersof twareprodukt zur Ermittlung 
von Bilddaten uber die Verteilung physikalisch und chemisch 
relevanter Daten in Untersuchungsobj ekten menschlicher oder 
10 tierischer Natur oder aus dem Bereich der Material- oder Si- 
cherheitsprufung. Insbesondere betrifft die Erfindung ein 
Auftrennen radiographischer Bilddaten in eine Verteilung der 
Materialdichte und eine Verteilung der Ordnungszahl im Unter- 
suchungsobj ekt . 

15 

Ergebnis aller radiographischen Verfahren, wie beispielsweise 
der Computer-Tomographie, der Mammographie, der Angiographie, 
der Rontgen-Inspektionstechnik oder vergleichbarer Verfahren, 
ist die Darstellung der Schwachung eines Rontgenstrahls ent- 

20 lang seines Weges von der Rontgenquelle zum Rontgendetektor . 
Diese Schwachung wird von den durchstrahlten Medien bzw. Ma- 
terialien entlang des Strahlengangs verursacht . Die Schwa- 
chung wird ublicherweise als der Logarithmus des Verhaltnis- 
ses der Intensitat der geschwachten zur Primarstrahlung defi- 

25 niert und bezogen auf ein Wegnormal als Schwachungskoef f i- 
zient des Materials bezeichnet. 

Eine Vielzahl radiographischer Untersuchungsgerate verwendet 
zur Darstellung der Schwachungsverteilung eines Rontgen- 

30 strahls in einem Untersuchungsgegenstand statt des Schwa- 
chungskoef fizienten einen auf den Schwachungskoef fizienten 
von Wasser normierten Wert, die GT-Zahl. Diese berechnet sich 
aus einem aktuell durch Messung ermittelten Schwachungskoef - 
fizienten ii und dem Referenz- Schwachungskoef fizienten Mh 2 o 

3 5 nach folgender Gleichung: 
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C = 1000 X [HU] (1) 

^H 2 0 

mit der CT-Zahl C in der Einheit Hounsfield [HU] • Fur Wasser 
ergibt sich ein Wert Ch 2 o=0 HU und furoLuft ein Wert C L = - 

5 1000 HU. 

Da beide Darstellungen ineinander trans formierbar bzw. acjui- 
valent sind, bezeichnet im folgenden der allgemein gewahlte 
Begriff Schwachungswert sowohl den Schwachungskoef fizienten m 

10 als auch den CT-Wert - Ferner werden im Sachzusainmenhang die- 
ser Erf indungsbeschreibung die Begriff e Material und Gewebe 
austauschbar verwendet. Es wird unterstellt, dass ein Materi- 
al im Kontext einer medizinisch angezeigten Untersuchung ein 
anatomisches Gewebe sein kann, und umgekehrt unter Gewebe in 

15 der Material- und Sicherheitspriifung ein beliebiges Material 
eines Untersuchungsobj ekts zu verstehen ist. 

Erhohte Schwachungswerte lassen sich entweder auf Materialmen 
hdherer Ordnungszahl, wie beispielsweise Calcium im Skelett 

2 0 oder Jod in einem Kontrastmittel zuriickfuhren, oder auf eine 
erhohte Weichteildichte, wie etwa bei einem Lung enkno ten . Der 
lokale Schwachungskoef fizient \x am Ort r ist abhangig von 
der in das Gewebe bzw. Material eingestrahlten Rontgenenergie 
E und der lokalen Gewebe- bzw. Materialdichte p entsprechend 

2 5 der folgenden Gleichung: 



fx = jut (E, r) = 



— (E) X p (r) (2) 



mit dem energie- und materialabhangigen Massenschwachungsko- 



30 effizienten 



(E) . 



Die energieabhangige Rontgenabsorption eines Materials, wie 
sie von seiner effektiven Ordnungszahl bestimmt wird, uberla- 
gert daher die von der Materialdichte beeinf lusste RSntgenab 
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sorption. Materialien bzw. Gewebe unterschiedlicher chemi- 
scher wie physikalischer Zusammensetzung konnen daher im 
Rontgenbild identische Schwachungswerte aufweisen. Umgekehrt 
kann dagegen aus dem Schwachungswert einer Rontgenauf nahme 
5 nicht auf die Materialzusamrnensetzung eines Untersuchungsob- 
jekts geschlossen werden . 

Eine korrekte Interpretation der somit eigentlich eher unan- 
schaulichen Verteilung der Schwachungswerte in einem mit ei- 

10 nem radiographischen Untersuchungsverfahren erstellten Ront- 
genbild kann im medizinischen Sektor meist nur aufgrund mor- 
phologischer Kriterien erfolgen und erfordert meist einen 
Radiologen mit jahrzehntelanger Erfahrung auf seinem Gebiet. 
Dennoch konnen in einigen Fallen Strukturen, die in der Bild- 

15 gebung einer Rontgenuntersuchung mit erhohten Schwachungswer- 
ten auf fallen, nicht klar klassif iziert werden. Bei spiel swei- 
se ist eine hilusnahe Verkalkung auf einer Thoraxiiber- 
sichts auf nahme nur schwer von einem orthograd zur Bildebene 
liegenden GefaS zu unterscheiden . Auch eine diffuse Kalkein- 

2 0 lagerung kann beispielsweise kaum von einer frischen Einblu- 

tung unterschieden werden. 

Auch in der Material- und Sicherheitspruf ung erganzt der Pru- 

> 

fer i.a. die Information der Darstellung einer Schwachungs- 
25 wert- Verteilung durch seine personliche Fachkenntnis und be- 
rufliche Erfahrung. Dennoch ist ihm z.B. ein sicheres Unter- 
scheiden einer kunststof f gebundenen Sprengstof fmischungen von 
einem nicht explosiven Kunststof f direkt aus einem Rontgenbild 
nicht moglich. 

30 

Hierfiir sind Verfahren zur Darstellung materialcharakteristi- 
scher Werte erf orderlich. W. Kalender et. al beschreiben in 
"Materialselektive Bildgebung und Dichtemessung mit der Zwei- 
Spektren-Methode f I. Grundlagen und Methodik, W. Kalender, W. 

3 5 Bautz, D. Felsenberg, C. SiiS und E . Klotz , Digit . Bilddiagn. 7, 

1987, 66-77 , Georg Thieme Verlag" ein Verfahren zur Basisma- 
terialzerlegung bei Rontgenauf nahmen . Das Verfahren basiert 
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auf dem Effekt, dass Materialien und Gewebe hoherer Ordnungs- 
zahl niederenergetische Rontgenstrahlung deutlich starker 
absorbieren als Materialien bzw. Gewebe niederer Ordnungs- 
zahl. Bei hoheren Rontgenstrahlenergien gleichen sich dagegen 
5 die Schwachungswerte an und sind vorwiegend eine Funktion der 
Materialdichte . 

Im Kontext dieser Beschreibung wird der Begriff Ordnungszahl, 
soweit nicht anders angegeben, nicht im strengen, elementbe- 
10 zogenen Sinn verwendet, sondern bezeichnet stattdessen eine 
effektive Ordnungszahl eines Gewebes , respektive Materials, 
die sich aus den ch.emisch.en Ordnungszahl en und Atomgewichten 
der am Aufbau des Gewebes bzw. Materials beteiligten Elemente 
berechnet . 

15 

Im von W. Kalender et . al vorgeschlagenen Verfahren werden 
die Rontgenschwachungswerte eines Untersuchungsobjekts mit 
Rontgenstrahlen niederer und hoherer Energie gemessen und die 
erhaltenen Werte mit den entsprechenden Ref erenzwerten zweier 

2 0 Basismaterialien wie beispielsweise Calcium (fur Skelettmate- 
rial) und Wasser (fur Weichteilgewebe) verglichen. Es wird 
angenommen, dass sich jeder Messwert als lineare Superpositi- 
on der Messwerte der beiden Basismaterialien darstellen 
lasst. Zum Beispiel kann fur jedes Element der bildlichen 

2 5 Darstellung des Untersuchungsobjekts aus dem Vergleich mit 
den Werten der Basismaterialien ein Skelettanteil und ein 
Weichgewebeanteil berechnet werden, so dass eine Transforma- 
tion der ursprtinglichen Aufnahmen in Darstellungen der beiden 
Basismaterialien Skelettmaterial und Weichteilgewebe resul- 

30 tiert. 

Die Basismaterialzerlegung bzw. das Zwei-Spektren-Verf ahren 
eignet sich damit zur Auftrennung bzw. Unterscheidung von 
anatomischen Strukturen in menschlichen und tierischen Gewe- 
35 ben mit stark unterschiedlicher Ordnungszahl. In der Materi- 
al- und Sicherheitsprufung konnte damit beispielsweise eine 
Auftrennung nach vordef inierten Materialarten, sogenannten 
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Materialklassen, erfolgen. Eine funktionelle Darstellung, die 
physikalische und chernische Charakteristiken der untersuchten 
Materialien bzw. Variationen dieser Charakteristiken inner- 
halb einer Materialart erkennen lasst, ist nicht Zielrichtung 
5 der Basismaterialzerlegung . 

Es ist daher Aufgabe der vorliegenden Erfindung ein Verfah- 
ren, eine Rontgenapparatur und ein Computersof twareprodukt 
vorzuschlagen, die auf der Basis radiograph! sch erstellter 
10 Absorptionsdaten eines Untersuchungsobjekts chemisch und phy- 
sikalisch relevante Information liber das Untersuchungsobjekt 
erstellen. 



Diese Aufgabe wird gelost durch ein Verfahren zur Ermittlung 

15 der Verteilungen von Dichte und Ordnungszahl in einem Unter- 
sucliungsobjekt mit Schritten zum Aufzeichnen einer ersten 
Verteilung von Rontgenabsorptionswerten des Untersuchungsob- 
jekts bei einem ersten Rontgenspektrum und zum Aufzeichnen 
zumindest einer zweiten Verteilung von Rontgenabsorptionswer- 

20 ten des Untersuchungsobjekts bei einem zweiten Rontgenstrahl- 
spektrum , zum Bestimmen einer ersten funktionalen Abhangig- 
keit eines ersten Rontgenabsorptionswertes der ersten Vertei- 
lung von Rontgenabsorptionswerten von Dichte und Ordnungszahl 
und zum Bestimmen einer zweiten funktionalen Abhangigkeit 

25 eines, dem ersten Rontgenabsorptionswert zugeordneten zweiten 
Rontgenabsorptionswertes der zweiten Verteilung von Rontgen- 
absorptionswerten von Dichte und Ordnungszahl, und 
zum Ermitteln eines Werts fur die Dichte und die Ordnungszahl 
aus einem Vergleich der ersten funktionalen Abhangigkeit mit 

30 der zweiten funktionalen Abhangigkeit und/oder weiteren funk- 
tionalen Abhangigkeit en . 



Der im Zusammenhang dieser Schrift verwendete Begriff 'Ront- 
genspektrum' besitzt eine weiter gefasste Bedeutung als nur 
3 5 die Spektralverteilung einer von der Rontgenquelle der Appa- 
ratur emittierten Rontgenstrahlung . Auf Seiten der Rontgende- 
tektoren werden unterschiedliche Spektralanteile einer. Strah- 
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lung mit unterschiedlichen Wirkungsgraden umgesetzt und somit 
verschieden gewichtet. Die daraus resultierende effektive 
Spektralverteilung wird in dieser Schrift als Rontgenspektrum 
bezeichnet . 

5 

Obige Aufgabe wird weiterhin gelost durch eine Rontgenappara- 
tur zur Ermittlung der Verteilungen von Dichte und Ordnungs- 
zahl in einem Untersuchungsobj ekt mit einer Rontgenquelle zur 
Emission von Rontgenstrahlung, einem Rontgendetektor zum 

10 Nachweis der von der Rontgenquelle emittierten Rontgenstrah- 
lung und zum Umsetzen der Rontgenstrahlung in elektrische 
Signale fur eine weitere Verarbeitung, und einer Signalverar- 
beitungseinrichtung fur das Verarbeiten der elektrischen Sig- 
nale des Rontgendetektors, wobei die Signalverarbeitungsein- 

15 richtung die Verteilungen von Dichte und Ordnungszahl in ei- 
nem Untersuchungsobj ekt auf Basis zumindest zweier, an der 
Rontgenapparatur mit unterschiedlichen Rontgenstrahlspektren 
aufgezeichneten Verteilungen einer Rontgenabsorption des Un- 
tersuchungsobj ekts nach einem erf indungsgemaSen Verfahren 

20 ermittelt. 



Ferner wird die Aufgabe gelost durch ein Computersof twarepro- 
dukt zur Ermittlung der Verteilungen von Dichte und Ordnungs- 
zahl in einem Untersuchungsobj ekt auf Basis von zumindest 
25 zwei bei unterschiedlichen Rontgenstrahlspektren aufgezeich- 
neten Bilddatensatzen der Verteilung von Rontgenabsorptions- 
werten in dem Untersuchungsobj ekt nach einem erf indungsgema- 
Sen Verfahren. 



3 0 Die vorliegende Erfindung ermoglicht die Berechnung der raum- 
liche Verteilung der mittleren Dichte p(F) und der effekti- 
ven Ordnungszahl Z(r ) aus einer Auswertung der spektral be- 
einflussten Messdaten einer Rontgenapparatur. Man erhalt 
hieruber neuartige Kontraste, insbesondere bezuglich der che- 

35 mischen und physikalischen Zusammensetzung des Untersuchungs- 
objekts. Diese, bislang Magnetresonanz-Systemen vorbehaltene 
funktionale Darstellung eines Untersuchungsobj ekts eroffnet 
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der Rontgendiagnostik wie der Rontgeninspektionstechnik eine 
Vielzahl neuer Anwendungen . 

Beispielsweise erlaubt die Darstellung der Verteilung der 
5 Ordnungszahl im Gewebe u.a. Einblicke in die biochemische 
Zusainmensetzung eines untersuchten Objekts, Kontraste auf - 
grund des chemischen Aufbaus in bisher dichtehomogen darge- 
stellten Organen, eine quantitative Bestimmung von Korperbe- 
standteilen wie z.B. Jod oder dergleichen und ein Herausseg- 
10 men tier en von Calcif izierungen basierend auf der Ordnungs- 
zahl. Die isolierte Dichtedarstellung eines Objekts gestattet 
eine genaue Schwerpunktbestimmung und Dichtevermessung von 
Objekten, wie sie u.a. beispielsweise bei Osteoporose vorge- 
nornmen werden • 

15 

Im Bereich der Sicherheitstechnik bedeutet dies eine zuver- 
lassigere Detektierbarkeit gefahrlicher Komponenten , insbe- 
sondere von Explosivstof fen. In der Materialpriifung eroffnet 
sich der Zugang zur quantitativen Untersuchung der Material - 
2 0 zusammensetzung und der Dichteverteilung in den Pruflingen. 

Weitere vorteilhafte Ausfiihrungen der vorliegenden Erfindung 
sind in den entsprechenden Unteranspruchen definiert. 

2 5 Vorteilhaft erfolgt die Bestimmung der funktionalen Abhangig- 

keit der Rontgenabsorptionswerte von Dichte und Ordnungszahl 
fur zumindest ein Rontgenstrahlspektrum mittels Referenzmes- 
sung an einer Eichprobe und/oder in Form einer Simulation auf 
der Basis eines physikalischen Modells, so dass eine anlagen- 

3 0 spezifische Abhangigkeit eines Rontgenabsorptionswertes von 

Dichte und Ordnungszahl erhalten wird. 

Vorzugsweise wird ein Umformen der Verteilungen der Rontgen- 
absorptionswerte in eine Verteilung der Dichte und eine Ver- 
teilung der Ordnungszahl fur jeden der zugeordneten Rontgen- 
3 5 absorptionswerte der ersten Verteilung der Rontgenabsorpti- 
onswerte und der weiteren Verteilungen der Rontgenabsorpti- 
onswerte auf der Grundlage der Ermittlung eines Wertepaares 
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fur Dichte unci Ordnungszahl so vorgenommen, dass das Werte- 
paar die bestimmten f unktionalen Abhangigkeiten der Rontgen- 
absorption von Dichte und Ordnungszahl fur das erste Rontgen- 
strahlspektrum und zumindest ein weiteres Rontgenstrahlspekt- 
5 rum erfiillt. Damit konnen Dichte und Ordnungszahl fur ein 

Bildelement einfach als Schnittmenge der funktionalen Abhan- 
gigkeiten der einander zugeordneten Rontgenabsorptionswerte 
der auf gezeichneten Verteilungen der Rontgenabsorptionswerte 
berechnet werden. 

10 

Vorteilhaf terweise weist das erste Rontgenspektrum eine Quan- 
tenenergie auf, die relativ zur Quantenenergie des zweiten 
Rontgenspektrums eine Rontgenabsorption durch den Photoeffekt 
begunstigt, so dass eine hohe Auflosung in der Bestiramung der 
15 Ordnungszahl en erhalten wird . 

In einer bevorzugten Ausfuhrungsf orm der vorliegenden Erf in- 
dung wird zum Verandern eines Rontgenstrahlspektrums fur das 
Aufzeichnen des Untersuchungsobj ekts eine Veranderung zumin- 

2 0 dest eines Betriebsparameters der Rontgenrohre vorgenommen, 

wobei die Rontgenquelle in einem ersten Betriebszustand ein 
erstes Ron tgens t r ahl spek t rum und in einem zweiten Betriebszu- 
stand ein davon verschiedenes zweites Ront gens trahl spek trum 
emittiert, so dass ein schneller Wechsel zwischen zwei Ront- 
25 genspektren ermoglicht ist. 

Ferner wird zum Verandern eines Rontgenstrahlspektrums fur 
das Aufzeichnen des Untersuchungsobj ekts vorteilhaf t eine 
Veranderung der Detektorcharakteristik vorgenommen, wobei der 

3 0 Rontgendetektor spektrale Teilbereiche der von der Rontgen- 

quelle empfangenen Ront gens trahlung in voneinander unabhangi- 
ge elektrische Signale umsetzt und hierbei ein simultanes 
Aufzeichnen von Verteilungen der Rontgenabsorption bei unter- 
schiedlichen Rontgenspektren zulasst . 

35 

Das Aufzeichnen zumindest einer Teilregion des Untersuchungs- 
obj ekts erfolgt ferner vorteilhaf t mit erhohter Rontgendosis 
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damit fur diese Teilregion eine zuf riedenstellende Auflosung 
der Dichte- und/oder Ordnungszahl -Verteilung erzielt wird 
ohne das gesamte Untersuchungsgebiet mit einer erhohten Roxit- 
gendosis zu belegen und z.B. einen Patienten einer zu hohen 
5 Strahlenbelastung auszusetzen. 



Zum Aufzeichnen einer erst en Verteilung einer Rontgenabsorp- 
tion des Untersuchungsobj ekts kann entsprechend einer vor- 
teilhaften Weiterentwicklung ein Material in den Strahlengang 

10 der Rontgenstrahlung eingebracht werden das zum Aufzeichnen 

einer zweiten Verteilung einer R5ntgenabsorption des Untersu- 
chungsobj ekts nicht in den Strahlengang der Rontgenstrahlung 
eingebracht wird. Damit kSnnen beim Aufzeichnen der ersten 
Schwa chungswertverteilung vor bzw. nach dem Untersuchungsob- 

15 jekt Spektralanteile der Rontgenstrahlung, die fur das je- 
weils verwendete Material spezifisch sind, herausgef iltert 
werden, so dass eine hohe Empf indlichkeit im Nachweis einer, 
diesem Material entsprechenden Substanz im Untersuchungsob- 
jekt erreicht wird. 

20 

Im Folgenden wird die vorliegende Erfindung anhand von Aus- 
fuhrungsbeispielen naher beschrieben, wobei auf folgende Fi- 
guren verwiesen wird, von denen 

25 Fig. 1 anhand einer Isoabsorptionslinie das Zustandekommen 

identischer Schwachungswerte /x bei Materialien un- 
terschiedlicher Zusammensetzung zeigt, 



Fig. 2 die Energieabhangigkeit der Rontgenschwachung fiir 
30 drei Elemente zeigt, 

Fig- 3a ein Funktionsschema der erf indungsgemaSen Berech- 

nungsmethoden zur Ermittlung von Isoabsorptionsli 
nien zeigt, 



35 
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Fig. 3b ein Ablauf diagramm der Transformation der Rontge- 

nabsorptionwerte in Werte der Materialdichte und 
Ordnungszahl zeigt , 

Fig. 4 zwei Isoabsorptionslinien einer Gewebeart bei zwei 

unterschiedlichen Ronfcgenstrahlspektren zeigt, 

Fig. 5 schematisch drei typische Rontgentechniken zeigt 

und 

Fig. 6 schematise!! verschiedene Methoden zur Variation der 

Rontgenstrahlspektren darstellt. 

Die Isoabsorptionslinie 11 des Diagramms 10 der Fig. 1 ver- 
bindet alle Wertepaare (p,Z) mit, bei einem definierten Ront- 
genstrahlspektrum identischem Schwachungswert \x bzw. C. Die 
Darstellung der Fig. 1 verdeutlicht , dass Information uber 
Art und Zusammensetzung eines Gewebes bzw. Materials nicht 
allein auf die Schwachungswerte eines Rontgenbildes gestutzt 
abgeleitet werden konnen . Ublicherweise benutzt ein Radiologe 
zum Identif izieren von Gewebearten im Rontgenbild seine Ana- 
tomiekenntnisse und sucht auf dieser Basis nach UnregelmaJSig- 
keiten. Zur Klarung der Identitat der UnregelmaEigkeiten ist 
ein Mediziner dann wiederum gezwungen auf Erf ahrungswerte und 
morphologische Kriterien zuruckgreif en . Ahnlich stiitzt sich 
ein Fachmann der Material- uns Sicherheitspriif ung zur Beur- 
teilung des radiographischen Befunds auf seinen beruflichen 
Erf ahrungsschatz . 

Rontgenstrahlung wird von unterschiedlichen Materialien und 
abhangig von der Energie der Rontgenstrahlung unterschiedlich 
stark geschwacht. Die Fig. 2 illustriert dies anhand der E- 
nergieabhangigkeit 2 0 des Massenschwachungskoef f izienten fur 
Wasser 23, Calcium 22 und Jod 21. Dies ist auf unterschied- 
lich wirkende Schwachungsmechanismen bei den verschiedenen 
Materialien zuriickzufuhren. Im diagnostisch relevanten Ener- 
giebereich der Rontgenstrahlung ist die Ron t gens chwachung im 
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wesentlichen auf die Absorption, verursacht durch den Photo- 
effekt und die, auf dem Comptonef f ekt beruhende Streuung zu- 
riickzufuhren . Die Absorption ist insbesondere bei niedriger 
Energie der Rontgenstrahlung und bei Geweben mit hoher Ord- 
5 nungszahl relevant. Die Streuung weist eine geringe Abhangig- 
keit von der Energie der Rontgenstrahlung auf und ist im we- 
sentlichen von der Elektronendichte, vermittelt uber die phy- 
sikalischen Dichte des Gewebes, abhangig. 



10 



15 



Die im Kontext dieser Beschreibung vereinfacht als Ordnungs- 
zahl titulierte ef fektive Ordnungszahl Z einer bestiminten 
Gewebeart errechnet sich aus den Ordnungszahl en Z± der am 
Aufbau beteiligten Elemente, deren Atomgewichte Ai und deren 
lokalen materialaquivalenten Dichten p± beispielsweise zu: 



Z = < 



v pi ^ 
L — piz ± 

i Ai 

pi 

£ — pi 

I i Ai 



(3) 



Fur reines Calcium erhalt man Z Ca =20 , fur Calciumhydrid ca. 
ZcaH2=16,04 und fur Wasser etwa Z H 2o=7,428. Die chemische oder 
20 auch biochemische Zusammensetzung eines Objekts kann daher 
sehr gut uber die Ordnungszahl Z erfasst we r den . 

Voraussetzung fiir eine Berechnung der Ordnungszahl- und Dich- 
teverteilung in einem Untersuchungsgebiet sind zumindest 

25 zwei, in der Aufnahmegeometrie identische, aber mit unter- 
schiedlicher Energie der angewandten Rontgenstrahlung er- 
stellte Rontgenaufnahmen des Gebiets. Bei Verwendung von mehr 
als zwei mit unterschiedlicher Rontgenstrahlenergie aufge- 
zeichneten Rontgenaufnahmen konnen die Z- und p-Auflosung 

30 verbessert werden, doch erhoht sich dadurch auch die Strah- 
lenbelastung. Im Falle der Untersuchung eines Patienten ist 
diese Moglichkeit daher nicht immer gegeben. 
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Ausgangspunkt der Umwandlung von Schwachungswert basierenden 
Bilddaten in Verteilungsbilder der Ordnungszahlen und der 
Material- bzw. Gewebedichte ist die Kenntnis der Isoabsorpti- 
onslinien fur jedes Rontgenspektrum einer Rontgenapparatur . 
5 Wie bereits erwahnt, ist hierbei unter Rontgenspektrum nicht 
der eng gefasste Begriff der Spektralverteilung einer von der 
Rontgenquelle der Apparatur emittierten Rontgenstrahlung zu 
verstehen, sondern ein erweiterter Begriff, der die unter- 
schiedliche Gewichtung unterschiedlicher Spektralbereiche des 
10 Emissionsspektrums der Rontgenrohre auf Seiten der Rontgende- 

tektoren beriicksichtigt . Ein gemessener Schwachungswert er- 

* 

gibt sich daher aus der direkten Schwachung des von der Ront- 
genrohre emittierten Strahlenspektrums und dern spektralen 
Wirkungsgrad des verwendeten Rontgendetektors . Beide Werte 
15 sind anlagenspezif ische Grofien und muss en entweder direkt 
oder indirekt mittels der Schwachungswerte von Eichproben 
ermittelt werden . Sie sind die Grundlage zur Berechnung der 
Isoabsorptionslinien. 

2 0 In der Fig. 3a sind drei Verfahren 3 00 zur Model lierung bzw. 
zur Berechnung einer Schar von Isoabsorptionslinien skiz- 
ziert, namlich eine theoretische Modellierung, eine experi- 
mentelle Bestimmung und eine theoretische Modellierung mit 
einer Kalibrierung der Kurven durch experiment ell bestimmte 

2 5 Parameter - 

Prinzipiell sind so viele Isoabsorptionslinien zu bestimmen, 
wie Schwachungswerte zum Abdecken der Spanne von Rontgen- 
schwachungen in den Rontgenauf nahmen erf order lich sind. Dabei 
30 ist nicht fur j eden theoretisch auftretenden Schwachungswert 
eine Isoabsorptionslinie zu berechnen; nicht errechnete Iso- 
absorptionslinien konnen bei Bedarf durch Interpolation oder 
andere geeignete Mi ttelungsverf ahren verfligbar gemacht wer- 
den. 

35 

Die Grundschritte der theoretischen Modellierung sind im lin- 
ken Ast des Ablauf schemas der Fig. 3a dargestellt. Im Schritt 
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S3 02 werden zunachst die Daten der fur eine Milage spezifi- 
schen Rontgenemissionsspektren S(E) mit den verfiigbaren Roh- 
renspannungen als Parameter eingelesen. Die Spektralvertei- 
lungen der Rontgenstrahlung konnen hierzu im Vorfeld experi- 
5 mentell fiir jede einzelne Rontgenanlage ausgemessen werden, 
oder es werden die fur einen speziellen Rontgenquellentyp 
charakteristischen Daten verwendet . Das Ermitteln der Detek- 
torapparatefunktion w(E) erfolgt in Schritt S303. Auch hierzu 
kann im Vorfeld eine genaue Vermes sung der Detektoranordnung 

10 vorgenommen werden oder aber es werden den Detektortyp cha- 
rakterisierenden Daten wie z.B. dessen spektrale technische 
Spezif ikation verwendet. Die Berechnung der Isoabsorptionsli- 
nien in Form von Kurvenscharen Ci(p,Z) bzw. Jii(p/Z) wird auf 
der Basis eines physikalischen Modells in Schritt S3 04 vorge- 

15 nommen, das fur jede relevante Kombination von S(E) und w(E) 
die Ron t gens chwachungen Ci bzw. fii fiir Materialien mit unter- 
schiedlichen Ordnungszahlen und bei unterschiedlichen Materi- 
aldichten nachbildet. 

2 0 Alternativ zur theoretischen Modellierung der Schritte S3 02 
bis S3 04 konnen die Kurvenscharen der Isoabsorptidnslinien 
auch experimentell ermittelt werden. Hierzu werden in Schritt 
S3 05 die Rontgenschwachungen von Eichmaterialien mit unter- 
schiedlicher Dichte und mittlerer Ordnungszahl in der Rfint- 

2 5 genapparatur bei verschiedenen relevanten Kombinationen von 

S(E) und w(E) gemessen. Die Messwerte bilden die Stiitzpunkte 
fiir die folgende Berechnung der Kurvenscharen von Isoabsorp- 
tionslinien Ci bzw. ix± in Schritt S306. 

Als weitere Alternative konnen die auf theoretischer Basis 

3 0 modellierten Kurvenscharen Ci bzw. Mi mit experimentell er- 

mittelten Rontgenschwachungswerten kalibriert werden. In 
Schritt S 307 werden die zum Eichen der theoretischen Kurven- 
scharen notwendigen Schwachungswerte wie oben fiir Schritt 
S3 05 beschrieben mit geeigneten Eichmaterialien bzw. Phanto- 
35 men in der Rontgenanlage gemessen. Im Unterschied zur rein 

theoretischen Modellierung der Schritte S3 02 bis S3 04 ist bei 
diesem Verfahren die exakte Kenntnis der Rontgenemissions- 
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spektren S(E) und w(E) nicht Voraussetzung sondern Parameter 
der theoretischen Modellierung der Kurvenscharen von Isoab- 
sorptionslinien Ci bzw. fj,± in Schritt S308. Das Kalibrieren 
der Kurven in Schritt S3 09 mit den in Schritt S3 07 experimen- 
5 tell ermittelten Eichwerten definiert schlieSlich Werte fur 
diese Parameter > die spezifisch fur die Ron t genemi s s i ons - 
spektren und Detektorapparatef unktionen der Rontgenapparatur 
sind* 

10 Mit der Ermittlung der Isoabsorptionslinien fur die erforder- 
lichen Rontgenschwachungswerte und Kombinationen von S(E) und 
w(E) sind die Voraussetzungen fur eine Transformation von 
Bilddaten, die Schwachungswerte der Rontgenstrahlung beim 
Durchgang durch ein Gewebe reprasentieren in Bilddaten, die 

15 eine Verteilung der Ordnungszahl bzw. der Materialdichte im 
entsprechenden Gewebe reprasentieren geschaffen. 

Je nach Auf gabenstellung konnen die drei Verfahren zur Isoab- 
sorptionslinienbestimmung auch gemischt verwendet werden . 

20 Beispielsweise konnen Werte, die experimentell nur ungenau 

oder nur mit grofiem Auf wand oder gar nicht zu ermitteln sind, 
mithilfe einer theoretischen Modellierung erganzt oder in 
ihrer Genauigkeit prazisiert werden. Die mit unterschiedli- 
chen Methoden erschlossenen Daten werden dann in Schritt S310 

25 zu einem einheitlichen Datensatz zusammengef asst und in 

Schritt S311 fiir die Bildtransformationen bereitgehalten. 



In der Fig. 3b ist das erf indungsgemaEe Transf ormationsver- 
fahren 32 0 dargestellt. Es stiitzt sich auf die nach einem der 
30 zuvor beschriebenen Verfahren 3 00 ermittelten und als Daten- 
satz in Schritt S321 bereitgehaltenen Kurvenscharen von Iso- 
absorptionslinien . 

Eine Transformation erfolgt bildelementweise . Im folgenden 
35 wird von einer Transformation einer Rontgenschwachungswert- 
verteilung basierend auf zwei bei unterschiedlichen Rontgen- 
strahlenergiespektren aber identischer Auf nahmegeometrie auf- 



WO 03/024331 



PCT/DE02/03055 



15 

genommen Ron tgenbi Idem ausgegangen. Dies ist die minimale 
Voraussetzung fur eine Durchfuhrung einer erfindungsgemafien 
Transformation. Jedoch konnen auch mehr als zwei Rontgenauf- 
nahmen bei mehr als zwei unterschiedlichen Energi evert eilun- 
5 gen der Ron t gens trahlung Verwendung finden. 

Die Auswahl eines zu trans formierenden Bildelements wird im 
Schritt S322 getroffen und im folgenden Schritt S323 werden 
die Schwachungswerte Ci bzw. [Xi fur dieses Bildelement aus 

10 dem ersten und C2 bzw. ix 2 aus dem zweiten Rontgenbild gele- 
sen. Im anschliefienden Schritt S324 erfolgt die Abfrage des 
flir die erste Rontgenaufnahme verwendeten Rontgenstrahlspekt- 
rums Si(E) und der Detektorapparatef unktionen wi(E) sowie der 
entsprechenden Werte S 2 (E) und w 2 (E) fiir das zweite Rontgen- 

15 bild. Diese Werte bilden die Parameter fiir eine nachfolgende 
Auswahl der den jeweiligen Schwachungswerten zuzuordnenden 
Isoabsorptionslinien. Die Spektralverteilungen S± (E) bzw. 
wi(E) konnen hierbei auch indirekt, z.B. uber eine Abfrage 
der verwendeten Rohrenspannungen Ui bzw. U 2 bzw. der Be- 

20 triebsparameter der Rontgendetektoren ermittelt werden. 

Im Schritt S325 werden aus dem in Schritt S321 bereitgehalte- 
nen Datensatz von Isoabsorptionslinien eine erste Kurve, wel- 
che die Bedingungen Ci bzw. 1x1 bei den Parametern Si(E) und 
25 wi(E) erfullt und eine zweite Kurve, welche die Bedingungen 
C2 bzw. jj.2 bei den Parametern S2 (E) und W2(E) erfullt ausge- 
wahlt. Ein Beispiel einer dergestalt erhaltenen ersten 11 und 
zweiten 41 Isoabsorptionslinie ist im Diagramm 40 der Fig. 4 
dargestellt . 

30 

Der Schnittpunkt 42 als Schnittmenge beider Kurven 11 und 41 
wird im Schritt S32 6 des Transf ormationsverf ahrens 32 0 be- 
rechnet. Der Kurvenschnitt 42 lasst sich z.B. durch eine lo- 
kale lineare Transformation oder mittels iterativer Schnitt- 
35 punktfindung ermitteln. Da die beiden Kurven 11 und 41 zwei 
unterschiedliche Schwachungswerte fur das selbe Bildelement 
und daher fur einen identischen Teilbereich eines untersuch- 
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ten Gewebes reprasentieren, imissen beide Schwachungswerte von 
der selben Material- bzw. Gewebeart verursacht sein. Die Ko- 
ordinaten (p,Z) des Kurvenschnittpunkts 42 geben daher die 
Materialdichte und die Ordnungszahl des dem Bildelement zuzu- 
5 ordnenden Gewebeteilbereichs wieder. 

SchlieSlich wird in Schritt S3 27 der so ermittelte Ordnungs- 
zahlwert Z in die Ordnungszahlverteilung als entsprechender 
Bildelementwert geschrieben, in Schritt S3 2 8 analog der er- 

10 mittelte Materialdichtewert p in die Dichteverteilung . Die 
Schritte S322 bis S328 werden fair alle verbleibenden Bild- 
punkte wiederholt, bis eine abschlieSende Bildausgabe in 
Schritt S3 2 9 erfolgen kann. Dabei kann der Schritt S3 2 4 iiber- 
sprungen werden, da die Spektralverteilungen S±(E) bzw. 

15 wi(E)fur alle Bildelemente eines Bildes identisch sind. 

Die Genauigkeit der Bestimmung der Dichte p und der Ordnungs- 
zahl Z wird im wesentlichen vom Rauschanteil der Schwachungs- 
werte und der spektralen Aufldsung der Messung beeinflusst. 
20 Der Rauschanteil kann z.B. uber die verwendete Rontgendosis 

beeinflusst werden. Die spektrale Auflosung ist eine Funktion 
des Abstands und der Breite der zur Messung verwendeten Ront- 
genspektren . 

25 Generell erhalt man bei Rontgenspektren relativ niedriger 

Energie eine Bevorzugung der Rontgenschwachung durch den Pho- 
toeffekt, bei Rontgenspektren relativ hoherer Energie eine 
Bevorzugung der Rontgenschwachung durch den Comptonef f ekt . 
Genauer ausgedruckt ist der Einfluss der Ordnungszahl auf die 

30 Rontgenschwachungswerte einer Aufnahme bei geringerer Rdnt- 

genenergie relativ grofier als bei hoherer Rontgenenergie . Der 
Einfluss einer Material- bzw. Gewebedichte auf die Rontgen- 
schwachungswerte verhalt sich dagegen genau umgekehrt . Vor- 
teilhaft wird daher zunachst ein erstes Rontgenspektrum so 

35 gewahlt, dass ein deutlicher Anteil an den ersten Rontgen- 

schwachungswerten vom Einfluss der Ordnungszahl en des unter- 
suchten Gewebes oder Materials herruhrt und ein zweites Ront- 
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genspektrum wird dann so gewahlt, dass die Dichten des Unter- 
suchungsobj ekts einen deutlichen Einfluss auf die zweiten 
Rontgenschwachungswerte nehmen . 

5 Fur die Computertomographie (CT) werden daher Quantenenergien 
der Ront gens t rahl spekt r en im Bereich zwischen 50keV und 
12 0keV bevor zugt , so dass ein ausreichender Energieabstand 
zwischen einem ersten und einem zweiten Rontgenspektrum ge- 
wahlt werden kann ohne die Rontgendosis in fur Patienten 
10 schadliche Bereiche erhohen zu mussen. In der Mammographie 
und der Mikro-CT konnen der untere Energiebereich sogar auf 
lOkeV bzw. IkeV festgelegt werden. 

Die Material und Sicherheitsprufung verwendet i.A. hohere 

15 Rontgendosen als bei der Untersuchung von Lebewesen ver tret- 
bar sind. Mit diesen hoheren Rontgendosen erhalt man ein 
meist wesentlich besseres Signal-zu-Rauschverhaltnis . AuSer- 
dem konnen auch wesentlich hohere Ron t gens t rahl energi en von 
bis zu lMeV verwendet werden. Damit kann, auch wenn die Ord- 

20 nungszahl bei einer bestimmten Rontgenstrahlenergie nur zu 
einem geringen Anteil zur Rontgenabsorption beitragt, ihr 
Wert dennoch mit guter Genauigkeit bestimmt werden, sofern 
bei der verwendeten Rontgendosis der Rauschanteil der Ront- 
genabsorptionswerte deutlich geringer ist als der relative 

25 Anteil der Ordnungs zahl an den entsprechenden Rontgenabsorp- 
tionswerten . Es kommt daher nicht so sehr auf die absolute 
energetische Lage der verwendeten Rontgenspektren an, sondern 
auf deren unterschiedliche energetische Lage zueinander. Die 
Wahl der absoluten Lagen wird dann von den verwendbaren Ront- 

30 gendosen und den nachzuweisenden Stoffen beeinflusst. 

Ausreichende Ordnungs zahl-Kontraste fur einen zuf riedenstel- 
lenden Einblick in die chemische Zusammensetzung eines Mate- 
rials oder die biochemische Zusammensetzung eines Gewebes 
35 werden uber die Rontgendosis und die verwendeten Rontgen- 
spektren eingestellt. Die notwendige Auflosung in der Dar- 
stellung der Ordnungs zahlvertei lung hangt von der diagnosti- 
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schen Auf gabenstellung ab. Fur eine quantitative Messung von 
typischen Korperbestandteilen wie z.B. Jod, Magnesium , Koh- 
lenstoff, Stickstoff, Sauerstoff, Natrium, Phosphor, Schwe- 
fel, Chlor, Kalium, Mangan, Eisen, Kobalt, Kupfer, Zink, Se- 
5 len, Molybdan oder dergleichen ist eine andere Z-Auflosung 
erforderlich als zum Erzielen eines chemischen Kontrasts ba- 
sierend auf z.B. dem Eisengehalt in einer Leber. Relativ ge- 
ring sind die Anforderungen an eine Z-Auflosung fur ein Her- 
aussegmentieren von Objekten aus einer Umgebung mit stark 
10 unterschiedlicher Ordnungszahl. Prinzipiell muss die Z- 

Auflosung besser sein als die Differenz der noch zu unter- 
scheidenden Ordnungszahl en. Gegenwartig ist eine Z-Auflosung 
von bis zu 0,1 erreichbar. 

15 Die Dichteverteilung kann zur Messung der Gewichtsverteilung 
eines Objekts gemafi der Gleichung 

m(r) = p(r) x Av (?) (4) 

20 ausgewertet und zur Schwerpunktsbestimmung verwendet werden. 
Anwendungen sind z.B. die Dichtevermessung von Knochen oder 
des Gehirns, aber auch die Bestimmung des tatsachlichen 
Schwerpunkts eines Prixf gegenstands . 

Zur Klassif izierung und Segment ierung von Korpermaterialien 
25 wie beispielsweise Leber, Nieren, wei£er und grauer Gehirn- 

masse etc. oder Materialklassen kann eine korabinierte Auswer- 
tung der Ordnungszahl- und Dichteverteilung vorgenommen wer- 
den, indem eine Fensterung der Dichte und der liber die Ord- 
nungszahl bestimmten Gewebe- bzw. Materialart vorgenommen 
30 wird. Zusatzlich zum Segmentieren eines z.B. malignen Gewebes 
kann mit dieser Methode beispielsweise gleichzeitig dessen 
Stickstof f gehalt bestimmt und fur die Dosisplanung in einer 
nachf olgenden Strahlentherapie verwendet werden. 

35 Um einen Vergleich mit der klassischen, Schwachungswerte dar- 
stellenden Radiographie zu ermoglichen, konnen die urspriing- 
lichen Bilder auch weiterhin als Standard- Rontgenbilder dar- 
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gestellt werden . Falls diese nicht mehr vorhanden sind, lasst 
sich aus der Ordnungszahl- und Dichteverteilung eines Unter- 
suchungsgegenstands riickwarts ein abwartskompatibles Standard 
Rontgenbild errechnen . 

5 

Bei dem im Vergleich. zur C ompu t er t omogr aphi e hohen Detektor- 
rauschen von Rontgenbildsystemen 1st die Dynamik der Rontgen- 
bildkontraste und somit die erzielbare Z- und p-Auflosung 
reduziert. Unter Umstanden wurde die Vorgabe einer hohen Z- 

10 bzw. p-Auflosung bei einer grofif lachige Untersuchung eine 

unzulassig hohe Dosisbelastung bedeuten. Nach den Ubersichts- 
messungen wird daher entsprechend den diagnostischen Erfor- 
dernissen in einem relativ kleinen Unter suchungsgebiet eine 
zweite Me s sung mit erhohter Dosis vorgenommen . In der Mairano- 

15 graphie empfiehlt sich diese Methode zur Untersuchung von 
Mikro-Calcif izierungen . 

In der Fig. 5 sind drei verschiedene Rontgensys t erne zur Il- 
lustration dreier Grundprinzipien von Rontgendurchleuchtungs- 
20 verfahren dargestellt. 

Fig. 5a zeigt ein Rdntgenbildsystem 500 zur projektiven Ab- 
bildung eines Untersuchungsobj ekts 503 auf einen zwei dimensi- 
onal en Flachendetektor 502. Ein Patient 503 oder ein Gegens- 
25 tand 503 befinden sich darin im allgemeinen wahrend einer 

Rontgenabbildung relativ zur Anordnun§ aus Rontgenrohre 501 
und Detektor 502 in Ruhestellung . 

Die Fig. 5b gibt als Beispiel einer Computertomographie- 
30 Vorrichtung schematisch das Prinzip eines Facherstrahlgerats 
5 00 mit rotierendem Detektorsystem 502 wieder. Die von der 
Rontgenrohre 501 facherformig emittierten Rontgenstrahlen 
durchleuchten das Untersuchungsobj ekt 503 aus einer bestimm- 
ten Winkelstellung und treffen schlieSlich auf eine zeilen- 
35 formige Anordnung diskreter Detektoren. Ein Auf nahmezyklus 
umfasst eine Vielzahl solcher Durchleuchtungen in unter- 
schiedlichen Winkelstellungen zum Unter suchungsgegens tand. 
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Bei dem in der Fig. 5c illustrierten Typ einer Rontgendurch- 
leuchtungsvorrichtung 500 wird der Untersuchungsgegenstand 
5 03 mit einer Transportvorrichtung wie z.B. einem Transport- 
5 band 506 durch den Strahlengang 504 der Anlage gef uhrt . Die- 
ser Typ von Rontgenapparatur wird hauptsachlich in der Si- 
cherheitsprufung verwandt . Als Detektor wird ein linienformi- 
ger, eindimensionaler , oft L-formiger Detektor verwendet . 
Beim Durchfahren durch den Strahlengang 5 04 wird der Gegens- 
10 tand 503 von den von der Rontgenrohre 501 emittierten Ront- 
genstrahlen zeilenweise abgetastet und aus den so erstellten 
Bildzeilen ein vollstandiges Durchleuchtungsbild des Gegens- 
tands zusammengesetzt . 

15 Die Fig, 6 zeigt verschiedene Moglichkeiten der Beeinf lussung 
des effektiven Rontgenstrahlspektrums, wie es sich aus dern 
Zusammenwirken der von der Rontgenrohre 501 emittierten 
S tr ahlungs ver tei lung 610 (611) und der spektralen Charakte- 
ristik der Detektorapparatefunktion 612 (620) ergibt. 

20 

Fig. 6a illustriert ein schnelles Umschalten der Rohrenspan- 
nung von einem Wert Ui auf einen Wert U 2 , so dass ein Wechsel 
zwischen zwei Rontgenemissionsspektren 610 und 611 erzielt 
wird- Vorteilhaft wird jede Teilaufnahme eines Aufnahmezyklus 

2 5 zuerst bei einer der beiden Rohrenspannungen und unmittelbar 

anschlieEend bei der anderen der beiden Rohrenspannungen auf - 
gezeichnet urn Artef akte durch Bewegungen z.B. eines Patienten 
wahrend des Aufnahmezyklus zu minimieren. Bei Rontgenappara- 
ten 500 mit kurzen Messzeiten wie z.B. Rontgenpro j ektionssys- 

3 0 temen kann auch jeder Aufnahmezyklus vollstandig bei jeweils 

einer der beiden Rohrenspannungen durchfahren werden. Statt 
einem Umschalten der Rohrenspannung konnen auch Filter oder 
Filtersysteme zwischen Rontgenrohre und Untersuchungsobjekt 
bzw. zwischen Untersuchungsobjekt und Rontgendetektor zur 
3 5 Beeinf lussung des Rontgenspektrums eingebracht werden. Der 

dargestellt Betrieb der Rontgenrohre bei einer Betriebsspan- 
nung Ui wird dann ersetzt durch einen Betrieb bei einer Stan- 
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dard-Betriebsspannung und Verwendung von keinem oder einem 
ersten Filter bzw. Filtersystem; ein Betrieb bei U2 wird 
durch Verwenden eines zweiten Filters bzw. Filtersys terns er- 
setzt . 

5 

Da eine Rontgenrohre 501 keine monoenergetische Strahlung 
sondern ein relativ breites Spektrum von Rontgenstrahlung 504 
emittiert, kann man, wie in der Fig. 6b gezeigt in einem 
Messvorgang mehrere Aufnahmen bei unterschiedlichen Rontgen- 

10 spektren erhalten indem man einen energiesensitiven Detektor 
620 verwendet . Dieser liefert fur i.a. nebeneinander angeord- 
nete Spektralbereiche getrennte Messsignale. Man erhalt 
Schwachungswerte somit simultan fur unterschiedliche, vonein- 
ander getrennte Spektralbereiche des Ront gen spekt rums also in 

15 einem Aufnahmezyklus eine, von der Aus fiih rungs form und Be- 

schaltung des Detektors definierte Anzahl von Rontgenbildern 
bei unterschiedlichen Strahlenergien . Solche Detektoren 620 
konnen als Schichtaufbau-Detektoren realisiert werden, wobei 
ausgenutzt wird, dass die Eindringtief e von Rontgenstrahlung 

20 in das Schichtsystem des Detektors von der Energie der Ront- 
gencjuanten bestimmt wird. Alternativ zu Schichtdetektoren 
konnen Quantenzahler als energiesensitive Detektoren verwen- 
det werden. 

25 Eine dritte Variante zur Beeinf lussung des Rontgenemissions- 
spektrums zeigt die Fig. 6c. Hierbei wird ein Material 630 
eines auf gabenspezif ischen Elements z.B. in Form schnell dre- 
hender Filterrader 630 zwischen Rontgenrohre 501 und Untersu- 
chungsobjekt 503 oder zwischen Untersuchungsobjekt 5 03 und 

30 Rontgendetektor 502 platziert. Jedes Element absorbiert auf 
charakteristische Weise Rontgenstrahlung in einem ebenfalls 
charakteristischen Spektralbereich. Befindet sich ein Fltigel 
des Filterrads im Strahlengang, so wird ein fur dieses Ele- 
ment typischer Spektralbereich im Strahlengang gedampft bzw. 

35 im guns tigs ten Fall eliminiert. Fur jede Teilaufnahme eines 
Aufnahmezyklus werden nun eine Messung ohne und eine Messung 
mit einem Filterradf lugel im Strahlengang vorgenommen . Ein 
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dem Filterradelement entsprechendes Element im Untersuchungs- 
objekfc kann dann nur in den Messungen ohne Filterradabdeckung 
des Strahlengangs einen Rontgenschwachungswert 63 2 aufweisen. 
Im anderen Fall ist der gemessene Schwachungswert 633 ideal - 
5 erweise gleich Null. Auf diese Weise kann mit hoher Prazision 
das Vorhanden- bzw. Ni ch tvorhandens ein eines bestimmten Ele- 
ments, wie beispielsweise Eisen oder Kupfer in einem Untersu- 
chungsobjekt nachgewiesen werden. Diese Nachweismethode ist 
insbesondere in der Mammographie und der Materialprufung von 
10 Bedeutung . 
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Patentanspruche 

1. Verfahren zur Ermittlung der Verteilungen von Dichte und 
Ordnungszahl in einem Untersuchungsobj ekt (503) mit Schritten 

5 zum Auf zeichnen einer ersten Verteilung von Rontgenabsorpti- 
onswerten des Untersuchungsobj ekt s (503) bei einem ersten 
Rontgenspektrum und zum Aufzeichnen zumindest einer zweiten 
Verteilung von Rontgenabsorptionswerten des Untersuchungsob- 
jekts (503) bei einem zweiten Rontgenstrahlspektrum , 

10 zum Bestimmen einer ersten funktionalen Abhangigkeit (11) 

eines ersten Rontgenabsorptionswertes der ersten Verteilung 
von Rontgenabsorptionswerten von Dichte und Ordnungszahl und 
zum Bestimmen zumindest einer zweiten funktionalen Abhangig- 
keit (41) eines, dem ersten Rontgenabsorptionswert zugeordne- 

15 ten zweiten Rontgenabsorptionswertes der zweiten Verteilung 

von Rontgenabsorptionswerten von Dichte und Ordnungszahl, und 
zum Ermitteln eines Werts fur die Dichte und die Ordnungszahl 
aus einem Vergleich der ersten funktionalen Abhangigkeit (11) 
mit der zweiten funktionalen Abhangigkeit (41) und/oder wei- 

20 teren funktionalen Abhangigkeit en. 

2. Verfahren nach Anspruch 1, 

dadurch gekennzeichnet 

dass die Bestimmung der funktionalen Abhangigkeit (11, 41) 
25 der Rontgenabsorptionswert e von Dichte und Ordnungszahl fur 
zumindest ein Rontgenstrahlspektrum mittels Ref erenzmessung 
an einer Eichprobe erfolgt. 

3 . Verfahren nach Anspruch 1 oder 2 , 
30 dadurch gekennzeichnet 

dass die Bestimmung der f unktionellen Abhangigkeit (11, 41) 
der Rontgenabsorptionswerte von Dichte und Ordnungszahl fur 
zumindest ein Rontgenstrahlspektrum in Form einer Simulation 
auf der Basis eines physikalischen Modells erfolgt. 

35 

4. Verfahren nach Anspruch 1, 2 oder 3, 
dadurch gekennzeichnet 
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dass ein Umformen der Verteilungen der Rontgenabsorptionswer- 
te in eine Verteilung der Dichte und eine Verteilung der Ord- 
nungszahl ftir jeden der zugeordneten Rontgenabsorptionswerte 
der ersten Verteilung der Rontgenabsorptionswerte und der 
weiteren Verteilungen der Rontgenabsorptionswerte auf der 
Grundlage der Ermittlung eines Wertepaares (42) ftir Dichte 
und Ordnungszahl so vorgenommen wird, dass das Wertepaar (42) 
die bestimmten funktionalen Abhangigkeiten der Rontgenabsorp- 
tion von Dichte und Ordnungszahl fur das erste Rontgenstrahl- 
spektrum und zumindest ein weiteres Ront gens trahlspekt rum 
erfullt • 

5. Verfahren nach einem der Anspruche 1 bis 4, 
dadurch gekennzeichnet 

dass das erste Rontgenspektrum eine Quant enenergie aufweist, 
die relativ zur Quant enenergie des zweiten Ron tgenspekt rums 
eine Rontgenabsorption durch den Photoeffekt begunstigt. 

6. Verfahren nach einem der Anspruche 1 bis 5, 
dadurch gekennzeichnet 

dass zum Verandern eines Ront gens trahlspekt rums fur das Auf- 
zeichnen des Untersuchungsobj ekts (503) eine Veranderung zu- 
mindest eines Betriebsparameters der Rontgenrohre und/oder 
eine Veranderung der Detektorcharakteristik vorgenommen wird. 

7. Verfahren nach einem der Anspruche 1 bis 6, 
dadurch gekennzeichnet 

dass das Aufzeichnen zumindest einer Teilregion des Untersu- 
chungsobj ekts (503) mit erhohter Rontgendosis erfolgt. 

8. Rontgenapparatur (500) zur Ermittlung der Verteilungen von 
Dichte und Ordnungszahl in einem Untersuchungsob j ekt mit 
einer Rontgenquelle (501) zur Emission von Rontgenstrahlung 
(504) , 

einem Rontgendetektor (502) zum Nachweis der von der Rontgen- 
quelle emittierten Rontgenstrahlung (504) und zum Umsetzen 
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der Rontgenstrahlung (504) in elektrische Signale fur eine 
weitere Verarbeitung, und 

einer Signalverarbeitungseinrichtung (505) fur das Verarbei- 
ten der elektrischen Signale des Rontgendetektors (502) 
5 dadurch gekennzeichnet, 

dass die Signalverarbeitungseinrichtung (505) die Verteilun- 
gen von Dichte und Ordnungs zah.1 in einem Untersuchungsobj ekt 
(503) auf Basis zumindest zweier, an der Rontgenapparatur 
(500) mit unterschiedlichen Rontgenstrahlspektren aufgezeich- 
10 neten Verteilungen einer Rontgenabsorption des Untersuchungs- 
obj ekts (503) nach einem der Verfahren der Anspruche 1 bis 7 
ermittelt . 



9 . Rontgenapparatur nach Anspruch 8 , 
15 dadurch gekennzeichnet, 

dass die Rontgenquelle (501) in einem ersten Betriebszustand 
ein erstes Rontgenstrahlspektrum (610) und in einem zweiten 
Betriebszustand ein davon verschiedenes zweites Rontgen- 
strahlspektrum (611) emittiert. 

20 

10. Rontgenapparatur nach Anspruch 8 oder 9, 
dadurch gekennzeichnet, 

dass der Rontgendetektor (502) spektrale Teilbereiche der von 
der Rontgenquelle (501) empfangenen Rontgenstrahlung (504) in 
25 voneinander unabhangige elektrische Signale umsetzt. 

11. Rontgenapparatur nach Anspruch 8, 9 oder 10, 
dadurch gekennzeichnet, 

dass zum Aufzeichnen einer ersten Verteilung einer Rontgenab- 
30 sorption des Untersuchungsobj ekts (503) ein Material (630) in 
den Strahlengang der Rontgenstrahlung (504) eingebracht wird 
und zum Aufzeichnen einer zweiten Verteilung einer Rontgenab- 
sorption des Untersuchungsobj ekts (503) kein entsprechendes 
Material in den Strahlengang der Rontgenstrahlung (504) ein- 
35 gebracht wird. 
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12 . Computersof twareprodukt zur Ermittlung der Verteilungen 
von Dichte und Ordnungszahl in einem Untersuchungsobj ekt 
(503) auf Basis von zumindest zwei bei unterschiedlichen 
Rontgenstrahlspektren auf gezeichneten Bilddatensatzen der 
5 Verteilung von Rontgenabsorptionswerten in dem Untersuchungs- 
objekt (503) nach einem der Verfahren 1 bis 7. 
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